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We propose an experimental scheme to probe the contribution of a single Dirac cone to the Hall 
conductivity as half-odd topological number sequence. In our scheme, the quantum anomalous Hall 
effect as in graphene is simulated with cold atoms trapped in an optical lattice and subjected to 
a laser-induced non-Abelian gauge field. By tuning the laser intensity to change the gauge flux, 
the energies of the four Dirac points in the first Brillouin zone are shifted with each other and the 
contribution of the single Dirac cone to the total atomic Hall conductivity is manifested. We also 
show such manifestation can be experimentally probed with atomic density profile measurements. 

PACS numbers: 37.10.Jk, 73.43.-f, 67.85.Lm, 71.10.Fd 



I. INTRODUCTION 



With easy production of graphene, this strictly two- 
dimensional material has evoked strong interest in con- 
dense matter and high energy physics in the past few 
years, for instance, the integer quantum Hall effect 
(QHE) has been observed in graphene [l[. Remarkably, 
the transverse Hall conductivity has an anomalous form: 
c^xy ~ 4(n+l/2)e^//i {n ~ 0,±1, ...), where n is the Lan- 
dau level (LL) index and the factor 4 accounts for double 
valley and double spin degeneracy. Note that the con- 
ductance sequence is shifted with respect to the standard 
QHE sequence by 1/2. This unusual sequence is now well 
understood as arising from the fact that the zero energy 
LL has half of the higher LLs degeneracy [2[ , and the Hall 
conductivity for each Dirac fermion is half-integer quan- 
tized as {n + l/2)e^/h, {n = 0,±1, ...)[!. Due to the 
presence of fermion doubling theorem [j|, which states 
that for a time reversal invariant system Dirac points 
must come in pairs, the contribution of the single Dirac 
cone to the total Hall conductivity of the graphene is hid- 
den. Even though one can find some lattice models which 
have an odd number of massless Dirac cones to avoid the 
fermion doubling theorem, the final Hall conductivity of 
the model is still integer quantized because of the hidden 
massive Dirac fermions [5|. In fact, the key reason lies 
in the famous Thouless-Kohmoto~Nightingale-den Nijs 
(TKNN) formula which guarantees that the total Hall 
conductivity is quantized as some integer numbers if the 
Fermi energy stays in the energy gap [6|. So far experi- 
ment is still lacking to demonstrate the half-integer Hall 
conductivity contributed from a single Dirac cone. 



Very recently, Watanabe and colleagues proposed a 
model to show that each Dirac cone does indeed con- 
tribute to the Hall conductivity with half-odd topologi- 
cal numbers [3] . In their work, they constructed a lattice 
model in graphene with complex second nearest-neighbor 
hopping integral among some particular lattice sites. The 
most interesting feature of the model is that the energies 
of the two Dirac points are shifted with each other, but 
their massless Dirac cones are still preserved. Through 
computing the topological number for the shifted Dirac 
cones, the half-integer contribution to the Hall conduc- 
tivity from single Dirac cone can be demonstrated. How- 
ever, as the authors themselves pointed out in their work, 
the lattice model is fictitious and it is extremely hard to 
realize the unusual second nearest-neighbor hopping pro- 
posed there. 

On the other hand, optical lattices populated with cold 
atoms offer a very promising alternative avenue to ex- 
plore rich fundamental phenomena of condensed matter 
physics [8|. For example, ultracold atoms in the optical 
lattices can be used to simulate the Bose-Hubbard model 
and the fractional quantum Hall effect, etc. [0]. In anal- 
ogy to the graphene, it is also proposed that the Dirac 
fermions can be simulated and detected by the ultracold 
atoms in the honeycomb optical lattice [10| and then 
demonstrated in details that the massive and massless 
Dirac particles can be realized in realistic conditions [Il| . 
In addition to the honeycomb lattices, the Dirac fermions 
can also be generated from a line-centered-square lattice 
or a lattice with T3 symmetry [12| . Furthermore, the 
Dirac cones can be realized with a square lattice sub- 
jected to an artificial non-Abelian gauge potential. The 
artificial gauge potentials acting on neutral atoms have 
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recently attracted considerable interest [13|. The gauge 
potential can lead to an effective magnetic field and pro- 
vide opportunities to simulate the physics of charges sub- 
jected to magnetic fields. Many recent works have shown 





FIG. 1: (Color online) The dispersion relation near the four 
Dirac cones in the FBZ for the different parameters a and 

/3. (a) a = /3 = 7r/2; (b) a = 7r/2 + 0.1, ^ = 7r/2; (c) 
a = 7r/2 + 0.05, ^ = 7r/2 - 0.05. 



that the Abclian and non-Abelian gauge fields can be re- 
alized using optical dressing |14| or laser-assisted hopping 
\W{ in optical lattice for observing the Hofstadter butter- 
fly [lal , atomic QHE [l7Hl9l , topological insulators 2G] , 
axion electrodynamics |2l| and Majorana fermion 22l |. 
Both Abelian and non-Abelian gauge fields in cold atoms 
have recently been generated experimentally [23| . 

Motivated by these ongoing developments, we propose 
an experimental scheme to manipulate single Dirac cone 
and probe its contribution to Hall conductivity as half- 
odd topological numbers with cold atoms in the opti- 
cal lattices subjected to an artificial non-Abelian gauge 
potential. We devise a distinct method to engineer the 
Dirac cone through changing the non-Abelian gauge fiux 
in optical lattice. In the 7r-fiux regime, there will be 
four zero energy Dirac points in the first Brillouin zone 
which leads to the anomalous Hall effect as in graphene. 
The key point in our paper is that, through changing the 
non-Abelian gauge fiux, the energies of the four Dirac 
points are shifted with each other, but their positions 
and massless Dirac cone characters remain intact. By 



manipulating the Dirac cones and calculating the corre- 
sponding topological Chern numbers, the contribution of 
the single Dirac cone to the total atomic Hall conductiv- 
ity as half-odd topological numbers is demonstrate indi- 
rectly. In our model, different gauge fluxes can be easily 
adjusted through tuning the laser field intensity. Further- 
more, we show that the atomic Hall conductivity could be 
detected with the standard density profile measurement 
used in cold atomic systems. 



II. SUMULATING AND ENGINEERING DIRAC 
CONE 

Let us consider a two-component ultracold fermion gas 
trapped in a two dimensional square optical lattice with 
the lattice sites (x = ma, y = na), where a is the lat- 
tice spacing, t is the nearest-neighbor hopping strength 
and the integers m and n are the position of the lat- 
tice site. We assume that the atomic collisions are ren- 
dered negligible through adjusting Feshbach resonance 
24| . The optical lattice setup is subjected to a gauge 
potential as A = {aay, Pa^ + 27r$m), where <Jx{y) is 
the Pauli operator of the two atomic components la- 
belled by two atomic sublevels of a hyperfine manifold, 
$ is the uniform Abelian magnetic fiux per unit cell and 
we set c = h = e = a^l. This gauge potential 
may be generated using a two dimensional optical su- 
perlattice based on laser-assisted spin-dependent hopping 
111 [H m HI. The gauge fluxes (a,/3) are determined 
by the Rabi frequencies of the two-photon off-resonant 
transition and can be tuned via changing the laser field 
intensity [11 [H [HI, [ll] . The tight-binding Hamiltonian 
of the system is written as 



J 



^ = -^ER+i 



,n^ •^m,n ^ ^m,n+l^ 



„t i(/3<T^-|-27r$ji 



H.c, 



(1) 



r 



where Cm,,n is a 2-component field operator defined 
on a lattice site {m,n). We firstly study the sys- 
tem in the non-Abelian gauge potential when the uni- 
form Abelian gauge flux $ = 0. By diagonalizing the 
above Hamiltonian, one can get the single- particle en- 



ergy spectra as Ek = hg (fc) ± Jh^ (k) + h'^ (fc), where 

ho{k) ~ ~2t{cos{a) cos{kx) + cos{f3) cos{ky)) , h^ (k) = 
— 2isin(/3) sin(fcy) a.nd hy (k) = — 2tsin(a) sin(fc2;). When 
the gauge fluxes a = /3 = tt/2 (vr-flux regime), this dis- 
persion relation leads to four Diracpoints around which 
the dispersion relation is linear [l^. The four Dirac 
points are located at ki = (0, 0), k2 = (0, tt), ka = (tt, 0), 
k4 = (tt, tt) in the first Brillouin zone (FBZ) with zero 
energies. Note that this 7r-fiux gauge potential is in the 



Abelian regime because its Wilson loop |VF| = 2 [T^. In- 
terestingly, the relative energies of the four Dirac points 
could be shifted through changing the gauge fiux from the 
Abelian to the non-Abelian regime {\W\ < 2). This can 
be seen clearly by expanding the Hamiltonian in Eq.(l) 
around the four Dirac points. Apart the conventional lin- 
ear dispersion relation related to the Dirac Hamiltonian, 
we find that the energies of Dirac points become 
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In the case of Abelian 7r-flux regime (a = /? = 7''/2), 
the energies of the four Dirac points are zero, which 
agrees with the previous description. However, if the 



gauge fluxes have been changed into the non-Abehan 
regime a = 7r/2 + Aa and (3 = 7r/2 + A/3, their ener- 
gies arc <5i?i(4) = =F2t(sin(AQ;) + sin(A/3)) and ^-£2(3) = 
=F2i(sin(Aa) — sin(A/3)) and would be shifted with re- 
spect to each other. As shown in Fig. 1, the energies 
of the four Dirac points vary with the gauge fluxes but 
the masslcss Dirac cone character is still kept. In partic- 
ular, we shift the energies of four and two Dirac points 
respectively in Fig. 1(b) and (c). 
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III. HALF-ODD TOPOLOGICAL NUMBER 



Now we turn to address the Hall conductivity in the 
Dirac cone regime. According to the famous TKNN ex- 
pression , the quantized value of the total Hall conduc- 
tivity is described as 
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where \ux) is the wave function of the energy band A and 
T^jgz i^ ^^"^ magnetic Brillouin zone (MBZ). The wave 
function \ux) forms a U(l) fibre bundle on the MBZ and 
the corresponding quantized value of the Hall conduc- 
tivity is the first Chern number which is a topological 
invariant of the U(l) bundle. If the Fermi energy Ep lies 
in one of the energy gaps, the total quantized value of 
the Hall conductivity is a sum of the first Chern num- 
bers corresponding to each energy band below the Fermi 
energy. In Fig. 2, using lattice gauge theory [2g|, we 
calculate numerically the Hall conductivity in the rcla- 
tivistic regime in the condition of Fig.l (same indexing). 
Here the laser-induced uniform abelian gauge flux plays 
the role of external magnetic field, which is used to break 
the time-reversal symmetry for getting the quantum Hall 
conductivity. In the presence of this magnetic field, the 
energy spectrum of the Dirac cones are splitted into many 
relativistic LLs. We also plot the relativistic LLs in each 
Dirac cone and the corresponding half-odd topological 
number sequences in the lower part of Fig. 2. 

In Fig. 2(a), the calculated result is for the original 
lattice model without shifting the Dirac cones. The Hall 
conductivity associated to each energy gap in the Dirac 
cone regime is { — 10, —6, —2, 2, 6, 10}, which has a step of 
4 as that in graphene. Compared with the usual Hall con- 
ductivity, it has an anomalous form a^y ~ 4(7i + 1/2) /h 
{n = 0,±1,..). This anomaly arises from the fact that 
each Dirac cone contributes to the total Hall conductiv- 
ity as a half-odd topological number sequence: n + 1/2. 
The factor 4 accounts for double valley and double spin 
degeneracy. As indicated in Fig. 2(a), the sum of the 
half-odd topological numbers over the four Dirac cones in 
each gap just corresponds to the calculated Chern num- 
ber. However, because the TKNN formula guarantees 
that the quantized values of the Hall conductivity for 
each gap are integer topological Chern numbers [y], the 
contribution of a single Dirac cone to the Hall conduc- 
tivity as half-odd topological numbers has been hidden. 



To manifest the contribution of the single Dirac cone, 
in Fig. 2(b-c), we shift the energies of the Dirac cones 
and expect that the effect of this shift on the Hall 
conductivity can give evidence of their own contribu- 
tion. In Fig. 2(b), the Hall conductivity is calcu- 
lated when one pair of Dirac cones is shifted equally 
toward positive energy and the other pair toward nega- 
tive energy, and the relativistic LLs in each shifted Dirac 
cone are also plotted in the lower part of the figure. 
The calculated result shows that the new Hall sequence 
is {-10, -8, -6, -4, -2, 0, 2, 4, 6, 8, 10}. By comparison, 
one also can find that the positions of the shifted rel- 
ativistic LLs coincide with the positions of the jumps 
of the new Hall plateaus, i.e. each Hall plateau corre- 
sponds to one of the new energy gaps induced by the 
above shift. This means that the appearance of the new 
Hall sequence is just because of the shifted Dirac cones. 
Furthermore, when the Fermi energy lies in one of the 
new gaps, as indicated in the lower part of Fig. 2(b), each 
shifted Dirac cone contributes to the total Hall conduc- 
tivity of the gap with a half-odd topological number. It is 
found that the sum of the four half-odd topological num- 
bers corresponding to each gap is equal to the numerical 
calculated Chern number of the gap. For example, for 
the central gap, the two positive (negative) shifted Dirac 
cones contribute equally to the Hall conductivity of the 
gap with a half-odd topological number —1/2 (1/2), their 
sum 2 X (—1/2 -I- 1/2) is just equal to the numerical cal- 
culated Chern number 0. As seen from this agreement, 
the half-odd contribution of single Dirac cone is mani- 
fested. In Fig. 2(c), we only shift one Dirac cone toward 
positive and one toward negative energy, and the other 
two ones keep intact. The induced new Hall sequence 
{-10, -9, -7, -6, -5, -3, -2, -1, 1, 2, 3, 5, 6, 7, 9, 10} be- 
comes more complicated. However, one can still find that 
the sum of the half-odd topological number sequences 
in the Dirac cones with and without shift agrees with 
the calculated Chern number sequence. In fact, differ- 
ent shifting cases make no difference on the underlying 
agreement. Note that this agreement is very surprising, 
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FIG. 2: (Color online) In the conditions of Fig. 1 (same index- 
ing), the Hall plateaus against the Fermi energy in the rel- 
ativistic regime is shown when <1> = 1/121. For comparison, 
in the lower part of each figure, we also plot the relativistic 
LLs and show the half-odd Chern number corresponding to 
each gap in the four Dirac cones respectively. Here we choose 
low magnetic flux ($ — 1/121) so that we can get more LLs 
in the Dirac cone for rendering their contribution to the Hall 
conductivity more pronounced. 



because there is no obvious reason to show that the su- 
perposition of field theory for the shifted Dirac cones can 
also give the same Chern numbers [3]. So, based on these 
agreements, we can demonstrate indirectly the contribu- 
tion of single Dirac cone as half-odd topological numbers. 
Due to TKNN, one can only have indirect signatures of 
half-odds, but our setup allows one to have such partial 
signs in a clearer and simpler fashion as compared to 
existing one Q- 



IV. DETECTION 

In the following, we explain how to detect the atomic 
Hall conductivity with the density profile measurement 
[271 [28j . The density profile of the trapped atoms can be 
measured through the time-of-flight (TOF) imaging with 
the light absorption |29[ instead of momentum distribu- 
tion, which is specific to two dimensional system [30j. 
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FIG. 3: (Color online) The atomic density profiles n{E''p)/A<b 
against the local Fermi energy for the Dirac Fermions in the 
condition of Fig. 1 (same indexing) and the solid (dash) line 
corresponding to <I> = 1/121 (2/121). 



Note that usually an external harmonic trap has been 
applied to confine the atoms in the optical lattice. If 
the trap potential varies slowly compared to the lattice 
spacing, the local density approximation (LDA) is well 



satisfied [31|. Under this approximation, the local Fermi 
energy is defined as Ep[r) = Ep — V{r), where Ep is the 
Fermi energy in the trap center and V{r) = rauP'r^ j^ is 
the global harmonic trap potential. The atomic density 
at temperature T is written as 



n{E'p) = I 



dk 



(2^) 



2/(fc,^f) 
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where /(fc,£'j,) = l/{cxp[(£'A; -£;j,)/T]-Hl} is the Fermi- 
Dirac distribution. If the local Fermi energy E^p (r) falls 
in a energy gap as plotted in Fig. 2, the atomic den- 
sity n{Ep) will depict a plateau. The Hall conductivity 
is related to the atomic density according to the Streda 
formula dxy = t^^Ie' t- -^^ can be seen from this for- 
mula, to measure the Hall conductivity, we need to choose 
two magnetic flux values, and identify the plateaus in 
both density profiles that correspond to the same gap. 
The density difference of the two plateaus divided by the 
magnetic flux difference gives the Hall conductance of the 
gap. In Fig. 3, at low temperature, we plot the density 
profiles of the atomic gas in the relativistic regime in the 
condition of Fig. 1 (same indexing) . Each case contains 
two density profiles for $i = 1/121 and $2 =2/121 and 
A$ = $2 — 'ft = 1/121. Through comparing the two 
density profiles for $1 and $2, the atomic Hall conduc- 
tivity can be derived as ha^y — (^2 — ni)/A$ [23, l28| . 



One can find that the derived atomic Hall conductivity 
from the calculated theoretical result in Fig. 3 agrees 
with the prediction in Fig. 2. 



V. DISCUSSION AND SUMMARY 

Finally, we give a brief discussion on the density profile 
measurement from an experimental point of view. Due 
to the external harmonic trap, the local Fermi energy 
decreases continuously from the trap center to edge. If 
we let the Fermi energy at the center of the trap higher 
than t, it automatically scans from t to —t. The local 
Fermi energy at the center of the trap can be increased 
by loading more atoms in the trap [S^l • It is not difficult 
to make the local Fermi energy higher than t, as it has 
been done in many previous experiment of fermions in 
optical lattices [32|- Moreover, there is no need to look 
for the exact value of the local Fermi energy as long as we 
do the comparison between plateau densities in the same 
energy gap. It means that, to a certain degree, the mea- 
surement result is insensitive to the influence of the ex- 
perimental noise on the parameters ($, a,/3), which also 
reflects the topological character of the quantized value 
of the Hall conductivity. Note that the measure result is 
also influenced by the temperature and the external trap 
frequency. The temperature is related to the visibility of 
plateaus in density profile. As shown in [23] , the temper- 



ature of the atoms should be cooled to the order of nk for 
observing discernible plateaus, which is easy than realiz- 
ing the degenerate Fermi gases. The external harmonic 
trap frequency determines the validity of the LDA. If the 
trap potential varies smoothly so that the confinement 
scale does not approach the scale of lattice constant, the 
LDA will be safe. Obviously, this requirement is within 
the current experimental technology. 

In summary, with cold atoms in a non-Abelian optical 
lattice, we have proposed an experimental scheme to sim- 
ulate the quantum anomalous Hall effect as in graphene. 
By tuning the laser intensity to change the gauge fluxes, 
we show that the contribution of a single Dirac cone to 
the Hall conductivity is a half-odd topological number 
sequence, and such topological numbers are detectable 
through the standard density profile measurement used 
in cold atom systems. 
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